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EFFECTS OF SELECTED TRAJECTORY PARAMETEFvS OH WEIGHT 


TRENDS IN THE SHUTTLE THERMAL PROTECTION SYSTEM 

By Donald M. Curry, James W. Tolin, Jr. , 
and Winston D. Goodrich 
Lyndon B. Johnson Space Center 


SUMMARY 


A computer program was developed to compute thermal protection sys- 
tem weights for given input entry conditions. An empirically derived 
heating model that includes laminar, transitional, and turbulent heating 
correlations as functions of angle of attack, free-stream velocity, and 
free-stream density was obtained from a first-order extrapolation of wind 
tunnel data to flight conditions. A unit area weight correlation for a 
reusable surface insulation was developed as a function of local heat 
load and entry time. The heating model, the thermal protection system 
correlation methodology, and a four-degree-of- freedom tra^ ctory inte- 
gration program were combined to develop a thermal protection system 
trajectory computer progre .. Analytical studies were performed for a 
typical delta wing Space Shuttle Orbiter to predict the thermal protec- 
tion system weight trends with entry parameter variations. The results 
of this specific study are described. 


INTRODUCTION 


The design and development of a minimum-weight thermal protection 
system (TPS) for the Space Shuttle are required to minimize weight and 
cost and to allow maximum payload capability. However, this development 
requires an optimum blend of material selection and system design, and 
this, in turn, is strongly dependent on the severity and duration of the 
aerothermodynamic environment experienced by the TPS. The procedure for 
obtaining au optimum blend of the design paramete ,i that minimize TPS 
weight within specific operational constraints is described in the fol- 
lowing paragraphs. Emphasis is placed on the gross trends of the re- 
sults as the operational constraints are allowed to vary, rather than 
on the absolute value of the resulting TPS weights. The results of this 
study should not be applied to absolute TPS weights for vehicles that do 
not conform to the specific design methodologies used in this study. 
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To blend the technologies required to perform this study, a computer 
program was developed to determine the optimum entry trajectory in rela- 
tion to the Space Shuttle TPS weight. The program can also be used to 
compute TPS weights for given input entry trajectories or for trajecto- 
ries defined by non-TPS criteria such as equilibrium glide, acceleration 
limit, or other prescribed control logic. 

To apply the TPS design and aerothermodynami c-environment predic- 
tion methodologies to the complete Shut tle/Orbi ter, the outer surface 
area of a typical delta wing Space Shuttle Orbit er (JSC-ObOA) was divided 
into 20 surface panels. Because of the assumed sy vietry of the vehicle 
and its environment, only one symmetrical half of the vehicle was neces- 
sary to perform this study. 

An average surface heat transfer rate environment was specified for 
each panel. Laminar, transitional, and turbulent heat transfer rates 
were specified as functions of vehicle angle of attack, free-stream 
velocity, and free-stream density. These heating rates were obtained 
from a first-order extrapolation of normalized wind tunnel data to flight 
conditions. Increases in local heating rates due to transitional and 
turbulent boundary layers were computed as functions of the Reynolds 
number per meter behind a normal shock. Critical values of the Reynolds 
number per meter behind a normal shock causing the onset of transitional 
and fully turbulent flow were assigned to each surface panel. This en- 
vironmental prediction procedure was sufficiently accurate for this 
study because the unit area TPS weight of a panel is proportional to the 
one-eighth power of the panel unit area heat load. 

Unit area weight correlations for the reusable surface insulation 
(RSI) were developed as a function of local heat load and entry time. 
Weights of the surface coating and adhesive bond system were included 
in the TPS weight calculation procedure. A TPS computer program was 
used to calculate the RSI thickness required to maintain prescribed 
bondline temperature limits after the initiation of the entry phase. 

This report describes the analytical heating model, TPS correla- 
tion methodology, and trajectory analyses that have been used in develop- 
ing the TPS trajectory computer program. 


SYMBOLS 


A 


C 


P 


surface panel area, meters 


specific heat, Joules/kilogram-kelvin 
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average slope of unit weight W/A curves (fig. 3) 

C£(t) the abscissa intercept (fig. 3) for W/A values of zero 

C^, values defined in equation (12) 

K thermal conductivity, watts/meter-kelvin 

2 

Q total heat, Joules/meter 

2 

q heat flux rate, watts/centimeter 

Re Reynolds number 

T temperature , kelvxn 

t heating time, seconds 

V velocity, meters/second 

2 

W/A weight of insulator per unit area, kilograms /meter 

a angle of attack, degrees 

Y flightpath angle, degrees 

c emissivity 

2 

U viscosity, newton-secona/meter 

3 

p density, kilograms/meter 

0 Stefan-Boltzmann constant 

Subscripts : 

1 location 

L local 

lam laminar 



k 


ns conditions behind normal shock 

ns/m per meter behind the normal shock 

ref reference, 0.305-meter-radius sphere 

s sea level 

trans transition 

turb turbulent 

00 free stream 


ANALYTICAL METHODS 


The TPS trajectory computer program contains aerothermodynamic heat- 
ing relationships and TPS weight correlations that have been incorporated 
into a four-degree-of-freedom trajectory computer program. The assump- 
tions used in these analyses and their validation are discussed in the 
following sections. 


Aerothermodynamic Heating Methodology 

The outer surface area of a typical delta wing Space Shuttle Orbiter 
configuration was diviued into 20 panels (fig. l). The areas and average 
aerodynamic heating data for each of the 20 panels are given in table I. 
The heating data given in table I represent laminae heating rates nor- 
malized by the stagnation point heating rate on a .305-meter-radius 
sphere. For the purpose of this study, these heating rates were obtained 
from test data that were averaged over the areas occupied by the panels. 
Transitional and fully turbulent boundary layer heating rates are pre- 
sented as a function of the Reynolds number behind a normal shock. This 
reference Reynolds number was used in an effort to account for the effects 
of real gas flow on the flow field at flight conditions. Critical values 
of the Reynolds number from the onset of transitional flow have been 
approximated from limited turbulent flow data on the Shuttle (ref. 1). 
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The heating rates and temperatures on any one panel will he deter- 
mined by the state of the boundary layer (i,e., laminar or turbulent), 
the angle of attack, and the free-stream density and velocity. Laminar 
heating rates are computed as 


^lam 



( 1 ) 


where 

4 ref 


values of 4 laB /4 ref 
from reference 2. 


are obtained from the data in -table I and 


4 


ref 


19 



( 2 ) 


When the Reynolds number per meter behind the normal shock ^ e ns / m ex ~ 

ceeds the value of (Re , ) on a given surface panel, the heating 

' ns ' m 'turb 

rate is assumed to be that from a turbulent boundary layer and is com- 
puted as 


^urb 


(Mf- h!>v 

\W\ 4 W 


( 3 ) 


where 4-) Brn /4 T . pf is obtained from table I, 
equation (2), and the approximation for 4 tur ^ ) /4 lam is represented by 


4 ref 


is computed with 
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2.5 


Re 


ns/m 


1 . 5(R e ns / m ) 


trans 


0.3 


.5) 


When the computed Reynolds number per meter behind the normal shock is 
between the values of (Re n ,) and /Re , \ boundary 

' ' 'trans ' 'turb 

layer is transitional; and the heating rates can be presented by a linear 
interpolation between the laminar and turbulent heating rates. 


*Hrans ~ ^lam + (\urb "* ^lam) 


Re ns/m ~ ( Re ns/m) 


trans 


Re 


as/x 'turb ' ^^KranB 


( 6 ) 


where 4 lam and q^.^^ 8X6 computed from equations (l) and (3). For 

example, the values of and (^ e ns / m ) from table I 

for panel 1 are 108 000 and 162 000, respectively. 

The trajectory computer program incorporates an option for select- 
ing one of two methods for computing the Reynolds number behind a normal 
shock. In one case, the Reynolds number is computed as 


Re 


ns/m 


p V p V 
ns ns _ 00 ” 


ns 


ns 


(T) 


where the viscojity u of the air is calculated behind the normal shock 
by using equilibrium real gas thermodynamic and transport properties. 

If knowledge of the flow properties behind the normal shock is not avail- 
able , the second program option is used to compute the Reynolds number 
for use in the heating rate calculations of equations (Is), (5), and (6). 
In this case, the Reynolds number per meter behind the normal shock is 
approximated by the relationship 


Re ns/m 


p V 

OO 00 


1.5 


= 172 


( 8 ) 
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Equilibrium panel surface temperatures can now be computed from the 
relationship 


T 


i 



(9) 


where q^ is the panel heating rate. Finally, panel heat loads 
are obtained from 


t* 


i 



dt 


( 10 ) 


Weight Correlation Methodology 


The TPS weights are dependent on the aerothermodynamic heati. 
vironment (laminar /turbulent flow) and the entry trajectory (orb-'t oj i.i- 
tude and inclination, entry range, and entry flightpath angle). To 
eliminate the necessity of performing detailed sizing of a TPS for each 
entry trajectory flown, a concise semiempirical method was developed to 
determine required RSI thickness over a wide range of entry trajectory 
variables. S. A. Mezines developed a similar technique (ref. 3) to es- 
timate ablator weights. For an ablative slab subjected to a co * 
heat rate, Mezines found that a correlation which used total ht J 

Q and heating time t of the entry trajectory could be used t„ . act 
ablation TPS weights. Since the Shuttle/Orbiter entry guidance operates 
in a constant heating rate mode during, the time of significant heating, 
the entry environment is similar to the constant heating test conditions 

3 1/8 

in Mezines' analysis. Thus, the same correlation parameter (Q^f 3 ) 

was used to develop a TPS-RSI unit weight correlation for Orbiter entries. 


Five Shuttle entry trajectories were used to develop a correlation 
cf RSI unit weight as a function of heat load and entry time. These 
trajectories represented entry ranges of 6802 to 15 hk2 kilometer.- , 
trajectory entry times of 31 to 56 minutes, and total reference heat 

Q ftp 

loads (0. 305-meter-radius sphere) of 7.65 * 10 to xU.02 x 10 J/m . 

The RSI unit weight correlations were computed -sing pressure and tem- 
perature dependent properties (taole II). For this study, the TPS con- 
figuration (fig- 2) consisted of P^I bonded to an aluminum structure and 
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was thermally modeled as a three-layer system (RSI* adhesive bond, and 
aluminum) . The model was assumed to be at 300 K at ohe start of entry 
with an adiabatic boundary condition imposed on the aluminum back surface 

3 1/d 

during entry. Curves of RST unit weight as a function of (Q^t^) 

were developed for the trajectories ummarized in table III. By using 
a TPS weight computation program (ref. 4), optimized RSI thickresses were 
calculated with a 442 K maximum aluminum structure temperature constraint 
and an RSI minimum material thickness of 0.648 centimeter. Mezines 1 

3 , 1/8 

correlation parameter (Q^t ) reduces the RSI unit weight data to a 

Ci 

series of almost straight lines from which an empirical relationship was 
developed. 

The basic relationship developed from these curves is 


(W/A) 


RSI 


■ c i(v 3 ) 1/a * 


At) 


(ii) 


where represents an average slope of the unit weight W/A curves 

plotted in figure 3 and ( t ) represents a relationship involving only 

the time from ent v v interface to the pov t at the end of entry when the 

2 

reference heating aoe is 22.7 kW/r.i . Assuming that 

C 2 (w . C 3 (t - t 0 ) * (12) 

where C^(t) is defined as the abscissa intercept for W/A values of 

zero, combining equations (11) and ( 12 ) yields the basic relationship 
for predicting RSI unit weight 
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Using the TPS model just described and the unit weight curves of fig- 
ure 3, the following relationship was obtained 

(W/A) RS i = 0.25633(^3) / - (t.036 x 10 _3 )(t - 1200) - 2U.U 

( 11 +) 

In computing W/A from equation (lh ) , the minimum W/A value on any 

2 

panel is limited to 1.528 kg/m (0.6U8-centimeter minimum RSI thickness). 

The total TPS weight for the Orbiter is determined by summing the 
weight of each surface panel as follows. 

20 

Total TPS weight A.(W/A) i (15) 

i=l 


In addition to the RSI unit weight, surface-coating and adhesive-bond 
unit weights must be included in the total TPS weight calculation of 
equation (15). 


Trajectory Analyses 

Trajectory model . - The data presented in the subsequent sections 
were obtained by using a four-degree-of-freedom trajectory integration 
program that incorporated a closed-loop entry guidance method. 

Down-range control is provided by modulating the roll angle , and 
cross-range control is provided by controlling the direction of the roll 
angle. The roll control logic is based on analytical range predictions 
and is divided into three phases: the constant heat r' ‘-.e phase, the 

equilibrium glide ranging phase, and the constant-g ranging phase. 

The constant heat rate phase begins at 0.1*9 m/sec^ and is used to 
stabilize the trajectory at pullout (y = 0), This phase is terminal! 
when the current velocity is less than 7^20 m/sec. The equilibrium 
glide ranging phase computes the roll angle (based on flying an equilib- 
rium glide profile) required to null the range errors. Control is trans- 
ferred to the constant-g ranging phase when the current load factor is 
equal to the const ant -load- factor level required to reach the target 
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(lU.7 m/sec ). The constant-g ranging phase computes the desired 
constant-g level required to null the range errors and converts the de- 
sired constant-g level into a roll angle command. Control is transferred 
to the transition logic at Mach 6. 

The data are based on the JSC-OUCA delta wing Space Shuttle Orbiter 
aerodynamics and a weight of 63 502 kilograms. Although the JSC-0U0A is 
no longer the baseline configuration, it is a representative Orbiter 
delta configuration. The 1962 standard atmosphere was used for all 
analyses. The trajectories were targeted to the NASA John F. Kennedy 
Space Center. 

Targeting .- The deorbit was nominally targeted to achieve the entry 
target line presented in figure k. The deorbit attitudes were selected 
to provide the deorbit velocity increment along the local horizontal. 

The nominal entry range was selected for each case by biasing the entry 
range from the minimum achievable entry range. The minimum achievable 
entry range is based on either a l6M K temperature limit at a point 
12.5 percent aft of the nose on the underside centerline (panel 1) or 
on a 2.5g load- factor limit, whichever results in the longer entry range. 
The constant angle-of-attack entries are constrained by the temperature 
limit , and the variable angle-of-attack entries are constrained by th** 
load- factor limit. 

The nominal target is biased so that the nominal entry range is 
1111 kilometers greater than the minimum achievable entry range. This 
bias was selected so that navigation, atmospheric, and aerodynamic dis- 
persions will not result in temperature or load-factor constraint viola- 
tions. The TPS weight, maximum load factor, and maximum panel 1 
temperature are presented in figure 5 as a function of entry range for 
a typical set of entry conditions. The nominal entry range is presented 
ir> figure 6 as a function of orbit altitude and inclination for constant 
and variable angle-of-attack entry profiles. In general, the trend is 
for range to increase as orbit altitude and inclination increase; this 
trend reverses, however, at lower inclinations. The reversal is caused 
by the oblateness effect associated with entry latitude at low 
inclinations 


RESULTS OF DELTA WING ORBITER STUDY 


Results of the JSC-0U0A delta wing Space Shuttle Orbiter study are 
presented to show the effects of entry trajectory parameters on the TPS 
weight. The TPS includes an RSI characterized by a 16 UU K surface tem- 
perature limit , except for the surface leading edges that have a higher 
temperature limit. 



11 


The TFS weight requirements are defined as a function of orbit in- 
clination, orbit altitude, entry range, and entry flightpath angle. Data 
are presented for the constant and variable angle-of-attack entry pro- 
files. These profiles were selected to minimize the TPS weight while 
providing a 2222-kilometer cross-range capability at 90° inclination. 

The TPS weights are based on the heating and TPS models discussed in the 
section entitled "Analytical Methods." 

The effects of orbit inclination, orbit altitude, entry range, and 
entry flightpath angle are presented in figures 7 to 11 for the constant 
angle-of-attack entry profiles. Data are presented for orbit altitudes 
that vary from 185 to 1296 kilometers and for inclinations that vary 
from 20.5° to l40°. Figures 7 and 8 present TPS weight change as a func- 
tion of orbit inclination and orbit altitude. The trend is generally 
as expected, with TPS weight increasing as inclination and orbit alti- 
tude increase; however, this trend reverses with lower inclination values. 
This is the same effect as that shown in figure 6. 

The reference TPS weight was obtained by using an entry trajectory 
that deorbited from an apogee of 500 kilometers at 55° inclination. As 
shown in figure 7, reducing the orbit altitude at this inclination to 
185 kilometers decreases the TPS weight by 352 kilograms , and increasing 
the orbit altitude to 926 kilometers increases the TPS weight by 
590 kilograms. Increasing the inclination from 55° to 90° for a 
500-kilometer orbit altitude increases the TPS weight by 3^0 kilograms. 

Figure 9 presents normalized contours of TPS weight in the orbit 
inclination plane. It must be noted that these weights are for the nom- 
inal range and nominal flightpath angle. The TPS weight must be increased 
to compensate for dispersions in these parameters; however, the trends 
noted are valid. 

Figure 10 illustrates the effect of entry range on the TPS weight 
for constant 30° angle-of-attack entry profiles. These data are for en- 
tries following deorbit from a 500-kilometer orbit altitude. The ref- 
erence point for the entry range shown in figure 10 is the nominal range 
presented in figure 6. For a 55° orbit inclination, the sensitivity 
of TPS weight to entry range is approximately 0.50 kg/km. 

Figure 11 presents the effect of variations in entry flightpath 
angle on the TPS weight For constant 30° argle-of-attack entry profiles. 
Data are presented for a 55° orbit inclination with 105- to 500-kilor.eter 
orbit altitudes. 

The TPS weight for variable angle-of-attack entry profiles is pre- 
sented as a function of orbit altitude, entry range, and entry flight- 
path angle in figures 12 to lU. Data are presented for inclinations 
of 55° and 90° and orbit altitudes of 185 to 926 kilometers. 
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Figure 12 shews the effect of inclination and orbit altitude on TPS 
veight. The reference zero weight point is for a 500-kilometer orbit at 
a 55° inclination. Entry at a 90° inclination for the same orbit alti- 
tude requires an additional 136 kilograms of TPS weight. Entry from a 
185-kilometer orbit Ititude requires 91 kilograms less TPS weight than 
entry from a 500-kij.ometer orbit at a 55° orbit inclination. Entry from 
a 926-kilometer orbit altitude requires 102 kilograms more TPS veight 
than entry from a 500-kilometer orbit altitude at this same inclination. 

Figure 13 presents the effects of entry range on TPS veight for 
55° and 90° orbit inclinations. The reference points for the entry 
range in these figures is the nominal entry range presented in fig- 
ure 6(b). The sensitivity of TPS veight to range is approximately 
0.71 kg/km for the variable angle-of-attack entry profile. Figure 1*+ 
illustrates the effect of entry flightpath angle on TPS veight for the 
variable angle-of-attack profile for a 55° orbit inclination and a 
500-kilometer orbit altitude. 

Table IV summarizes the veight sensitivity dett -mined by this trade- 
off study. 


CONCLUSIONS 


A computation procedure has been developed to combine an empirically 
derived heating model, thermal protection system (TPS) correlation meth- 
odology, and trajectory analysis. This program can be used to compute 
TPS weights for given input entry trajectories or for trajectories de- 
fined by non-TPS criteria such as equilibrium glide, acceleration limit, 
or a prescribed control logic. 

ibis computation procedure was applied to the Lyndon B. Johnson 
Space Center 0U0A delta wing Orbiter to show the effects of entry tra- 
jectory parameters on the TPS veight trends, and TPS weight sensitivity 
coefficients were obtained. 

Although the computation procedure developed cannot be used to pre- 
dict exact TPS weights, gros3 trends of the TPS veight can be predicted 
as vehicle operational constraints are allowed to vary. 


Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, January 18 , 197** 
986-15-31-0L-72 
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TABLE I.- AERODYNAMIC HEATING RATE DATA FOR THE OSC-OLOA DELTA WING ORBITER PANELS 


Panel number 

(a) 

. 2 

Area, m 

Angle of 
attack, deg 

^lam/^ref 

(Be , ) 

\ ns/mJ 
' trans 

(b) 

1 

2U.8 

30 

0.221 

108 000 



53 

.260 


2 

53. L 

30 

.091 

66 GOO 



53 

.170 


3 

T5 .2 

30 

, 

.065 

39 000 



53 

.098 


k 

81 

30 

.065 

20 000 



53 

.085 


5 

Lo.l 

30 

.051 

108 000 



53 

.017 


6 

57.6 

30 

.031 

66 000 



53 

.OIL 


T 

38.? 

30 

.020 

13L 000 



53 

.008 


8 

**1.3 

30 

.005 

13L 000 



53 

.005 


9 

26.9 

30 

.026 

108 000 



53 

.026 


10 

6L.2 

30 

.008 

108 000 



53 

.003 


11 

88.6 

30 

.009 

108 OOO 



53 

.009 


12 

103.5 

30 

.00L 

108 000 



53 

.OOL 


13 

60.3 

30 

.020 

13L 000 



53 

.023 


1L 

*♦9.5 

30 

. 10L 

108 000 



53 

.156 


15 

48.7 

30 

.065 

67 000 



53 

.085 


16 

51.5 

30 

.033 

13L OOO 



53 

.033 


IT 

1*0.9 

30 

.020 

13L 000 



53 

.020 


18 

UL.9 

30 

.021 

ljL OOO 



53 

.021 


19 

27 

30 

.025 

108 000 



53 

.016 


20 

2L 

30 

.020 

108 000 



53 

.016 

1 


a See figure I. 

b Ncte: (Re . \ ® 1.5 * (Re . ^ 

' n9/ra Wb * na/m 'tr* 
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TABLE II.- THERMAL PHYSICAL PROPERTIES 
Constant values: p ■ 2bO kg/m 3 ; e * 0,85] 

(a) Thermal conductivity 



Thermal conductivity K, Vl/m-K, at a pressure of 

K 

1*337 * 10 1 S/m 2 

1.337 - 10 2 B/a 2 

1.337 « 10 3 H/m 2 

1.337 * 10 U H/a 2 

1.013 * 10 5 B/m 2 

0 

1.7068 * 10' 4 

1.8992 x 10“** 

2.8007 * IQ* 4 

3.2093 * 10* 4 

3.51 » io‘ u 

U22 

2.1035 

2.3078 

3.1012 

3.^17 

3.726 

6kh 

2.70U5 

2.6968 

3.5098 

U.6878 

5.229 

1033 

U.50T5 

L.6878 

6.01 

9.015 

10.00 

1367 

6.7312 

6.9115 

10.217 

IU.06 

15.63 

16UU 

9.1352 

9.3756 

13.^62 

18.87 

20.91 


(b) Specific heat 


Temperature, K 

Specific heat C , 
J/kg-K p 

200 

300 

300 

632 

hOO 

829 

500 

980 

5bh 

1101 

700 

1172 

1000 

1231 

1U00 

1323 

1605 

13U0 
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TABLE III.- THE JSC-OUOA ORBITER ENTRY TRAJECTORIES 


Case 

Trajectory- 
time , 
min: sec 

Down range, 
km 

Maximum 

6 

q « maximum, 
^ref ’ 

W/m 2 

Q ref total » 
J/m 2 

H 

ra 

9 692 

1.5 

9.56 x 10 5 

O 

10.80 x 10 



6 9U7 

1.5 

15.39 

8.62 


mEM 

7 875 

1.5 

11.70 

9.45 



15 442 

1.5 

7.20 

14.02 

1 5 

KB 

6 802 

2.2 

17.09 

7.65 


TABLE IV.- WEIGHT SENSITIVITY OF THE TPS 


Parameter 

Constant 
angle of attack 

Variable 
angle of attack 

Nominal mission range, km . - . 

9112 

6852 

TPS weight sensitivity to entry 
range , kg/km 

0.50 

0.71 

TPS weight sensitivity to orbit 



altitude kg/km 

1.1 

0.24 

TPS weight sensitivity to 



Q 

inclination, kg/deg 

14.7 

3.8 

TPS weight sensitivity to entry 



8L 

flightpath angle, kg/deg . . . 

912 

1263 


^or 55° inclination, 500-km orbit altitude. 

Above 500-km orbit altitude for 55° inclination. 
C Above 55° inclination for 500-km orbit altitude. 
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Figure 1.- Surface panels of a JSOOUOA delta wing Orbiter. 
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Figure 2.- Schematic of the TPS configuration. 
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Figure 3.- The TPS unit weight as a function of total heat load and time. 
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Figure 5»- Thermal protection system weight and maximum load factor 
compared to range for typical entry (inclination, 55°; variable 
angle of attach; orbit altitude, 500 kilometers). 
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Figure 7.- 
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Figure 8.- Change in TPS weight as a function of deorbit apogee altitude 
(angle of attack, 30°; nominal range). 
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